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ABSTRACT 

A  simple  analytical  relationship  between  a  ship- 
hull  form  and  its  steady  far-field  Kelvin  wake  is  obtained 
by  considering  the  low-Froude-number  limit  of  the 
Neumann-Kelvin  theory.  In  particular,  this  relationship 
predicts  the  occurrence  of  a  sharp  peak  in  the  amplitude 
of  the  waves  in  the  far-field  Kelvin  wake  at  an  angle,  a, 
from  the  ship  track  that  is  smaller  than  the  Kelvin-cusp 
angle  of  1 9°  1/2  for  a  hull  form  which  has  a  small  region 
of  flare  and  is  wall  sided  elsewhere,  if  the  Froude 
number  is  sufficiently  small.  An  explicit  relationship 
between  the  angle,  cp,  between  the  ship  track  and  the 
tangent  to  the  ship  mean  waterline  in  the  region  of  flare 
and  the  corresponding  “wave-peak"  angle  a  in  the 
Kelvin  wake  is  obtained.  For  instance,  this  relationship 
predicts  the  occurrence  of  a  sharp  peak  in  wave 
amplitude  at  an  angle  a  in  the  Kelvin  wake  equal  to  14° 
for  a  hull  having  a  small  region  of  flare  within  which  the 
waterline-tangent  angle  cp  is  approximately  equal  to  either 
30°  or  74°.  This  theoretical  result  may  explain  the  bright 
returns  that  have  sometimes  been  observed  in  SAR 
images  of  ship  wakes  at  angles  smaller  than  the  Kelvin- 
cusp  angle.  The  low-Froude-number  asymptotic  analysis 
of  the  Neumann-Kelvin  theory  presented  in  this  study 
also  predicts  that  the  wave-resistance  coefficient  is  0(F2), 
where  F  is  the  Froude  number,  for  a  ship  form  with  a 
region  of  flare.  OfF4)  for  a  ship  form  that  is  wall  sided 
everywhere  and  has  either  a  bow  or  a  stern  (or  both)  that 
is  neither  cusped  nor  round,  and  OfF6)  for  a  wall-sided 
ship  form  with  both  bow  and  stern  that  are  either  cusped 
or  round. 


1.  INTRODUCTION 

It  has  been  observed,  see  for  instance  Fu  and  Holt 
(1982)  and  McDonough  et.  al.  (198$),  that  SAR 
(Synthetic  Aperture  Radar)  images  of  ship  wakes 
sometimes  reveal  bright  returns  along  rays  at  angles  from 
the  track  of  the  ship  smaller  than  the  Kelvin-cusp  angle 
of  19°l/2.  A  plausible  explanation  for  these  surprising 
observations  was  proposed  by  Scragg  (1983)  who 
considered  a  simple  ship  bow  form  with  large  flare  for 
which  he  found  that  the  zeroth-order  slender-ship 
approximation  to  the  far-field  wave-amplitude  function 
given  in  Noblesse  (1983)  predicted  a  sharp  peak  in  the 
value  of  the  amplitude  of  the  divergent  waves  at  an  angle 
from  the  track  of  'he  ship  equal  to  approximately  half 
the  bow  entrance  angle.  This  numerical  result  of  Scragg 
was  confirmed  by  Barnell  and  Noblesse  (1986)  who  also 
found  that  the  peak  in  the  amplitude  of  the  divergent 
waves  becomes  sharper  as  the  value  of  the  Froude 
number  decreases,  and  thus  suggested  that  the  occurrence 
of  a  sharp  peak  in  the  amplitude  of  the  far-field  Kelvin 
waves  was  a  large-flare  low-Froude-number  feature. 

The  numerical  studies  of  Scragg  and  of  Barnell 
and  Noblesse  are  based  on  two  simple  approximations  to 
the  far-field  wave-amplitude  function,  namely  the  Michell 
thin-ship  approximation  for  which  no  peak  was  found 
and  the  zeroth-order  slender-ship  approximation  which 
exhibited  a  peak  as  was  already  noted,  so  that  it  is  not 
dear  from  these  studies  whether  a  more  realistic 
mathematical  model  for  the  far-field  wave-amplitude 
function,  such  as  that  provided  by  the  Neumann-Kelvin 
theory,  would  also  predict  the  occurrence  of  peaks  in  the 


amplitude  of  the  far-fteld  Kelvin  waves.  Furthermore,  the 
numerical  results  obtained  by  Scragg  and  by  Barnell  and 
Noblesse  correspond  to  a  particular  ship  form  and  thus 
provide  little  physical  insight  into  the  origin  of  the 
predicted  peak  in  the  amplitude  of  the  fa.-field  Kelvin 
waves,  specificallv  the  manner  in  which  such  a  peak  is 
related  to  the  shape  of  the  ship  hull. 

A  complementary  analytical  study  of  the  low- 
Froude-number  limit  of  the  Neumann-Kelvin  theory  for 
an  arbitrary  ship  form  is  thus  presented  here.  This 
asymptotic  analysis  of  the  Neumann-Kelvin  theory 
provides  a  simple  analytical  relationship  between  a  ship- 
hull  form  and  its  steady  far-field  Kelvin  wake.  In 
particular,  this  relationship  predicts  the  occurrence  of  a 
sharp  peak  in  the  amplitude  of  the  waves  in  the  far-fteld 
Kelvin  wake  at  an  angle,  a,  from  the  ship  track  that  is 
smaller  than  the  Kelvin-cusp  angle  of  19°  1/2  for  a  hull 
form  which  has  a  small  region  of  flare  and  is  wall  sided 
elsewhere,  if  the  value  of  the  Froude  number  is 
sufficiently  small.  A  simple  explicit  relationship  between 
the  angle,  <p,  between  the  ship  track  and  the  tangent  to 
the  ship  mean  waterline  in  the  region  of  flare  and  the 
corresponding  wave-peak  angle  a  in  the  Kelvin  wake  is 
given  and  depicted  in  Figure  3b.  For  instance,  this  figure 
predicts  the  occurrence  of  a  sharp  peak  in  the  amplitude 
of  the  divergent  or  transverse  waves  at  an  angle  o  in  the 
Kelvin  wake  equal  to  14°  for  a  hull  having  a  small  region 
of  flare  within  which  the  waterline-tangent  angle  cp  is 
approximately  equal  to  30°  or  74°,  respectively.  This 
analytical  result  may  explain  the  bright  returns  that  have 
sometimes  been  observed  in  SAR  images  of  ship  wakes 
at  angles  smaller  than  the  Kelvin-cusp  angle. 

The  low-Froude-number  asymptotic  analysis  of  the 
Neumann-Kelvin  theory  presented  in  this  study  also 
predicts  that  the  wave-resistance  coefficient  is  ()(F2), 
where  F  is  the  Froude  number,  for  a  ship  form  with  a 
region  of  flare,  0(F4)  for  a  ship  form  (hat  is  wall  sided 
everywhere  and  has  either  a  bow  or  a  stern  (or  both)  that 
are  neither  cusped  nor  round,  and  G(Fft)  for  a  wall-sided 
ship  form  with  both  bow  and  stern  that  are  either  cusped 
or  round 

More  precisely,  the  low-Froude-number  asymptotic- 
approximation  (50)  to  the  far-field  wave-amplitude 
function  Ktt)  shows  that  the  main  contributions  to  the 
function  Ktt)  stem  from  several  particular  points  on  the 


mean  waterline.  These  are  the  bow  and  the  stern,  on  one 
hand,  and  (usually  but  not  always)  one  (or  several) 
point(s)  of  stationary  phase.  Indeed,  the  number  of  these 
points  of  stationary  phase,  and  their  position  on  the 
waterline,  depend  on  the  value  of  t  and  on  the  shape  of 
the  waterline.  The  first  two  terms  in  the  low-Froude- 
number  asymptotic  expansions  for  the  contributions  KB  s 
of  the  bow  and  stern  and  the  contributions  Kt  of  the 
points  of  stationary  phase  in  equation  (50)  are  given  by 
equations  (51)-(55)  and  (J9)-(62),  respectively.  The 
second-order  terms  in  these  asymptotic  expansions  are 
defined  by  complex  expressions.  However,  the  first-order 
terms  provide  simple  approximations  defined  explicitly  in 
terms  of  the  geometrical  characteristics  of  the  hull  and 
the  velocity  components  in  the  tangential  directions  T and 
TT  x  T  to  the  hull  (see  Figure  2).  In  particular,  these  low- 
Froude-number  asymptotic  expansions  show  that  the 
contributions  KB  and  Ks  of  the  bow  and  stern  are  0(1) 
except  if  the  bow  or  stern  is  cusped  or  round,  in  which 
case  we  have  KB  s  =  0(F2).  The  contribution  of  a  given 
point  of  stationary  phase  is  0(l/F),  and  thus  is 
dominant,  if  the  hull  has  flare  at  this  point;  otherwise, 
that  is  if  the  hull  is  wall  sided  at  the  point  of  stationary 
phase,  its  contribution  is  O(F).  The  low-Froude-number 
approximation  (50)  also  shows  that  we  have  K(t)  = 
0(l/t5)  as  t  —  °°  In  fact,  this  result  is  valid  for  any 
value  of  the  Froude  number. 

The  latter  result  implies  that  the  lines  along  which 
the  steepness  of  the  short  divergent  waves  in  the  far-field 
Kelvin  wake  takes  given  large  values,  say  1/7  and  1/15, 
are  parallel  to  the  ship  track,  as  was  found  in  Figure  21 
of  Barnell  and  Noblesse  (1986)  by  using  the  Michell  thin 
ship  approximation  for  a  simple  ship  form.  The 
Neumann-Kelvin  theory  therefore  predicts  that  the  far- 
field  Kelvin  wake  contains  three  distinct  regions:  (i)  a 
narrow  constant-width  inner  region  bordering  the  track 
of  the  ship  where  no  divergent  gravity  waves  can  exist, 

(ii )  an  outer  region  where  the  usual  transverse  and 
divergent  waves  are  present,  and  (iii)  an  intermediate 
region  at  the  boundary  between  the  inner  and  outer 
regions  where  short  steep  divergent  waves  can  be  found 
In  reality,  surface  tension  and  possibly  also  viscosity 
must  evidently  be  taken  into  account  in  the  vicinity  of 
the  track  of  the  ship 


2.  STATIONARY-PHASE  APPROXIMATION  TO 
THE  FAR-FIELD  KELVIN  WAKE 

The  far-fleld  Kelvin  wake  may  be  conveniently 
analyzed  in  terms  of  the  nondimensional  far -field 
coordinates  (4.rj.4>  =  (X.Y.Z)  g/U2,  where  g  is  the 
gravitational  acceleration.  U  is  the  speed  of  advance  of 
the  ship,  and  (X,Y,Z)  are  dimensional  coordinates.  The 
mean  free  surface  is  taken  as  the  plane  {  =  0,  with  the 
{  axis  pointing  upwards,  and  the  4  axis  is  chosen  along 
the  track  of  the  ship,  that  is  in  the  ship  centerplane,  and 
pointing  towards  the  bow.  The  origin  of  the  system  of 
coordinates  is  placed  within  the  ship. 

Equations  (4),  (3a),  (7)  and  (8)  in  Barnell  and 
N<  lesse  (1986)  yield  the  following  expression  for  the 
elevation  (  of  the  free  surface  at  a  sufficiently  large 
distance  behind  the  ship,  such  that  nonlinearities  may  be 


neglected: 

niU.n)  =  ReJ^°(E++E_)K(t)<l+t2)',2dt,  (1) 

where  K(t)  is  the  far-field  wave-amplitude  function  and 
E„  is  the  exponential  function  defined  as 

E±(t;4,o)  =  exp  (i(6±(t;o)],  (2) 

with  the  phase-function  and  the  parameter  o  defined 
as 

0t(t;o)  =  (l*ot)((  +  t2)1'2,  (3) 

o  =  r,/<-{).  (4) 


For  a  ship  with  port-  and  starboard-symmetry,  as  is 
considered  here,  the  Kelvin  wake  is  symmetric  about  the 
ship  track  n  =  0.  We  may  then  restrict  the  analysis  of 
the  Kelvin  wake  to  the  quadrant  rj  >  0,  4  <  0  and 
assume  a  >  0.  Let  a  be  the  angle  between  the  track  of 
the  ship  and  the  line  joining  the  origin  of  the  system  of 
coordinates  to  the  observation  point  (4.tj).  We  thus  have 

a  =  tan  'o.  (5) 

A  far-field  asymptotic  approximation,  valid  as  4 
-<*>,  to  the  wave  integral  (1)  may  be  obtained  by 
applying  the  method  of  stationary  phase,  as  is  well 
known.  The  result  of  this  classical  asymptotic  analysis 
may  be  found  in  Barnell  and  Noblesse  (1986),  for 
instance.  Specifically,  equations  (28),  (24a),  (20b),  (25a-d) 
and  (26a-d)  in  this  reference  yield 

<-4),/J«4o)  ~  Re 


with  C<  s  <  tan"  ,(2“v2)  -  !9°28\ 

where  t±(a)  are  the  values  of  t  for  which  the  phase- 
function  0  +  (t;o)  is  stationary,  ©^(o)  =  8  +  (t±;o) 
represents  the  corresponding  values  of  the  phase-function 
0  +  (t;o),  and  A ±(a)  is  the  function  defined  as  A±(o)  = 
(2/*)l/2(l-(.tt2)1/2/|4:e';(t:1:;<*))l/2.  The  functions  tt(o), 


©t(a)  and  A x(a)  are  given  by 

t±(a)  =  [1±(1-8o2)1/21/4o,  (7) 

0t(<»)  =  (3t<1-8o2)1/2) 

|l  +4o2±(1  -8o2)1/2]i/2/2i/28o,  (8) 

At(o)  =  (I  +4o2±(I -8o2)|/2|3/42|/4/ 

4n‘  V/2(l -8o2),/4,  (9) 

where  we  have 

o  =  tana  (10) 

as  is  given  by  equation  (5). 


The  far-field  asymptotic  approximation  (6)  shows 
that  the  wave  pattern  at  any  point  (4,a),  with  4  «  -  1 
and  0  <  a  <  19“28 ' ,  consists  in  two  elementary  plane 
progressive  waves,  so-called  transverse  and  divergent 
waves,  as  is  most  well  known.  The  wavelengths  At  and 
the  directions  of  propagation  pz,  measured  from  the 
track  of  the  ship,  of  the  transverse  (A  ,/)_ )  and 
divergent  (A,  ,/)„)  waves  at  an  angle  a  from  the  track  of 
the  ship  are  given  by 

At(o)  -  2'  216no2/|3*(l  So2)1  :) 

[1  -  4o2±(I  -  8c2)1,2)1  2.  (II) 

pja)  =  cos  1  {21  22o/|I  +4o2±tl  -8o2)'  2)‘  2}.  (12) 

We  have 

I  >  A  /2n  >  2/3  A  ,  /2n  ^  0  and  (13) 

0  <  P  _  <  sin  1  ( 1  / 3 1  2>  St  35°  16  <  ft  ,  <  90°  (14) 

Equation  (6)  shows  that  the  amplitudes  of  the 
transverse  and  divergent  waves  in  the  Kelvin  wake  are 
asymptotically  equal  to  A,  K(tt)|/(  -4)1  2  as  4  —  ~°° 

The  steepnesses,  say  s,.,  of  these  waves  then  are  given  by 
st  =  At|K(tt)|/At(  -  4),/2.  We  then  have 

(  4)1  '2sJ4 -a)  ^  S1|K(tr)1  as  4  —  (15) 

where  the  functions  S,(o)  are  defined  as  S l(o)  = 
AI(o)/A1(a).  Equations  (9)  and  (II)  then  yield 

S.(o>  =  |3+( I  -  802)1  2)  |l  -  4o2±( I  -  802)1'2]1'2 


(A  K(t  )  exp  [i(4©  -  "24)) 

-t-  AtK(t  J  exp  |i(4©,  +n/4)|}  as  4  - 


|)  +  4o2±( I  802)1  2|3  4-'2'  464rt'  V  2(l  -  802)1  4.  (16) 


(6) 


Equations  (7),  (8),  (9),  (10)  and  (16)  show  that  we 


have 

t_  -  0,  8.  -  I,  A.  -  (2/«)l/2 

(17a,b,c) 

and  S.  -  l/23/V/2  for  a  -  0, 

(17d) 

r  +  rv  \/2a,  8+~  l/4< r,  A  +  ~  l/2*1/2a3/2 

(18a, b,c> 

and  S  +  ~  1/16h3/2o7/2  as  a  -»  0, 

(18d) 

t±  =  1/21/2  and  8±  =  3(3/2)l/2/4 

(19a, b) 

fora  =  tan“'(l/23/2)  19°28', 

A±  ~  3/n,/26l/4(l  -8c2)174  and 

(19c) 

SI  ~  9/4n3/261/4(l-8o2)|/4 

(19d) 

as  a  —  tan  '(l/23/2). 

The  stationary-phase  values  t±(o),  the  phase-functions 
6 t(a),  the  amplitude-functions  A t(a)  and  the  steepness- 
functions  SI(a)  are  depicted  in  Figures  la  and  lb,  which 
correspond  to  the  transverse  and  divergent  waves. 


a 


respectively.  Fig.  la  -  The  Stationary-Phase  Value  t  _  (a),  the  Phase- 

,  .  ,  .  ,  Function  8  (»),  the  Amplitude-Function  A  (o)  and 

It  may  be  seen  from  these  figures  and  from  the  Steepness-Function  S .  (a)  Corresponding  to  the 

equations  (19c)  and  (19d)  that  the  amplitude-functions  Transverse  Waves  in  the  Kelvin  Wake 


A±(a)  and  the  steepness-functions  S±(ai  are  singular  at 
the  Kelvin  cusp  line  o2  =  1/8 ,  a  =  tan  '  '(l/23/2)  ~ 
19°28\  in  accordance  with  the  well-known  fact  that  the 
asymptotic  approximation  (6)  is  not  uniformly  valid  at 
the  boundary  of  the  Kelvin  wake.  A  complementary 
asymptotic  approximation,  expressed  in  terms  of  Airy 
functions,  valid  at  and  near  the  Kelvin  cusp  line  is  given 
in  Ursell  (1960)  for  the  particular  case  of  the  Kelvin  wave 
pattern  due  to  a  concentrated  pressure  point  at  the  free 
surface.  However,  Ursell’s  more  complex  asymptotic 
approximation  *  ’)  not  be  considered  here  because  the 
simple  asymptotic  approximation  (6)  is  little  affected  by 
the  weak  singularity  (1  8o2)  1  4  for  points  ({.«)  inside 

the  Kelvin  wake  and  not  too  near  the  cusp  line  a  cm 
19°28',  in  which  we  are  mostly  interested  in  this  study. 

Equation  (II)  yields 

A  t  ~  8m»2  as  a  -*  0.  (20) 


This  approximation  anu  the  approximation  (18b)  show 
that  there  are  an  infinite  number  of  divergent  waves  with 
indefinitely  shorier  wavelength  in  the  vicinity  of  the  track 
of  the  ship,  as  is  well  known.  It  may  be  seen  from 
Figure  lb  and  equations  (I8a,c,d)  that  the  stationary- 
phase  value  t^(o),  the  amplitude-function  A  +  (o)  and  the 
steepness-function  St(o)  are  unbounded  in  the  limit  a  -» 
0.  Equations  (I8a,c,d)  yield 


a - ► 

Fig.  Ib  -  The  Stationary-Phase  Value  t^fer),  the  Phase- 
Function  6  +  (a),  the  Amplitude-Function  A  +  (a)  and 
the  Steepness-Function  S + (<r)  Corresponding  to  the 
Divergent  Waves  in  the  Kelvin  Wake 


A+  ~  (2,'ir),/2t  + 3/2  and  (2 Ja) 

S+  ~  t  +  7/2/2,/2tt3/2  as  t+  -»  <«>.  (21b) 


m 


Equations  (21a)  and  (6)  show  that  the  amplitude  of  the 
divergent  waves,  given  by  A  +  |K(t  +  )|/(-t),/2.  vanishes 
at  the  track  of  the  ship  if  we  have 

t3/2K(t)  -*  0  as  t  -»  <*>.  (22) 

Furthermore,  it  is  shown  in  Barnell  and  Noblesse  (1986) 
that  the  asymptotic  approximation  (6)  is  uniformly  valid 
at  the  track  of  the  ship  if  condition  (22)  is  verified.  It 
may  be  seen  from  equations  (21b)  and  (15)  that  the 
steepness  of  the  divergent  waves,  given  by 
S+  |K(t  +  )|/(  -4)|/2,  is  unbounded  at  the  track  of  the 
ship  if  we  have 

t7/2|K(t)|  —  °°  as  t  —  (23) 

Both  conditions  (22)  and  (23)  can  be  satisfied 
simultaneously  if  we  have 

|K(t)|  ~  l/tM  as  t  -»  00  with  3/2  <  gi  <  7/2.  (24) 

In  summary,  the  asymptotic  approximation  (6) 
expresses  the  far-field  wave  pattern  of  a  ship  at  a  point 
(4.o),  with  (  «  -  I  and  0  <  a  <  19°28\  as  the  sum  of  a 
transverse  wave  and  a  divergent  wave.  The  phase 

of  these  two  waves  are  defined  explicitly  in 
terms  of  (  and  a.  and  their  amplitudes 
AI(cr)K(t;t)/(  -{)'  2  are  given  by  the  product  of  the 
functions  A±(o)/(  -()'  2,  which  are  also  defined 
expi  citly  in  terms  of  (  and  o,  and  the  far-field  wave- 
amplitude  function  K[t,(o)l,  which  depends  on  the  speed 
(Froude  number)  and  the  shape  of  the  ship.  The  far-field 
wave-amplitude  function  i  now  considered. 


3  BASIC  FXPRESSIONS  FOR  THE  FAR  FIELD 
WAVE  AMPLITUDE  FUNCTION 

The  far-field  wave-amplitude  function  K(t)  may  be 
conveniently  defined  in  terms  of  the  nondimensional 
near  field  coordinates  (x.y.e)  =  (X,Y,Z)/L,  where  L  is 
the  length  of  the  ship  In  the  Neumann  Kelvin  theory, 
the  function  K(t)  is  given  b-  the  sum  of  an  integral 
around  the  mean  waterline  of  the  ship  and  an  integral 
ovc  ne  mean  wetted  hull  surface  Specifically,  for  a 
ship  with  port  and  starboard  symmetry,  as  is  considered 
here,  the  function  K(tl  may  be  expressed  in  the  form 

K<>)  =  K  .  (t)  *  K  (t),  (25) 

where  the  functions  K,  are  given  by 

F2K.(t)  =  ^  f'.'V  *  nA*j  *  iv;p4M,  dl 

<■  v:  j cxp(v2p2/)l  ,  (nx  +  v:p‘4n  . )  da.  (26) 


as  is  shown  in  Noblesse  (1983).  In  this  equation,  F  is  the 
Froude  number  defined  as 

F  =  U/(gL)12,  (27) 

v  is  its  inverse,  that  is 

v  =  1/F,  (28) 

p  is  defined  as 

p  =  (I +t2)l/2,  (29) 

and  Ex  represents  the  exponential  function 

E*  =  exp  (  -  iv2p(xj:ty).  (30) 

Furthermore,  c  and  h  represent  the  positive  halves  of  the 
mean  waterline  and  of  the  mean  wetted-hull  surface, 
respectively.  The  unit  vector  tangent  to  c  and  pointing 
towards  the  bow  is  denoted  by~t(tx,ty,0),  and  n)nx,ny,ny) 
is  the  unit  vector  normal  to  h  and  pointing  into  the 
water,  as  is  indicated  in  Figure  2  The  term  n±  is  defined 


n±  =  ny  +  i(nxttn>,)/p. 


Figure  2  -  Definition  Sketch  for  a  Single-Hull  Ship  with 
Port  and  Starboard  Symmetry 


Also,  dl  and  da  repiesent  the  differential  elements  of  arc 
length  of  c  and  of  area  of  h,  respectively.  Finally. 

$  =  4(x)  represents  the  nondimensional  disturbance 
potential  $  =  <t>  UI  at  the  integration  point  Ton  c  or  h. 
♦  represents  the  derivative  of  4  in  the  direction  of  the 
tangent  vector  fto  c.  and  4j  is  the  derivative  of  4  in  the 
direction  of  the  vector  if «  C,  which  is  tangent  to  h  and 
pointing  downwards  as  is  shown  in  Figure  2. 

Equations  (25)  and  (26)  express  the  far-field  wave 
amplitude  function  K(t)  in  terms  of  the  value  of  thr 
Froude  number  and  the  form  of  the  mean  wetted-hull 
surface,  as  was  noted  previously.  More  precisely, 
equation  (26)  expresses  the  function  K„(t)  as  the  sum  of 
a  line  integral  around  the  mean  waterline  and  a  surface 
integral  over  the  mean  wetted-hull  surface  of  the  ship 
Furthermore,  these  integrals  involve  the  disturbance 
velocity  potential  4  in  t heir  integrands  The  relationship 
between  t ho  far  field  wave  paticin  of  a  ship  and  its  form 


VVe’ ""V-  T- . 


and  speed  thus  is  a  complex  one,  which  offers  little 
physical  insight.  However,  analytical  approximations  for 
the  waterline  and  the  hull  integrals  in  equation  (26)  can 
be  obtained  in  the  limiting  case  when  the  value  of  the 
Froude  number  is  sufficiently  small.  This  low-Froude- 
number  asymptotic  approximation  for  the  function  K(t) 
is  now  obtained 

4.  LOW-FROUDE  NUMBER  LAPLACE 

APPROXIMATION  TO  THE  FAR-FIELD  WAVE- 
AMPLITUDE  FUNCTION  IN  TERMS  OF  A 
WATERLINE  INTEGRAL 

For  large  values  of  v  =  1/F,  or  more  generally  of 
vp,  the  exponential  function  exp(v2p2z)  in  the  hull- 
surface  integral  in  equation  (26)  vanishes  rapidly  for 
negative  values  of  z.  Therefore,  only  the  upper  part  of 
the  mean-wetted  hull  surface  h  yields  a  significant 
contribution  in  the  low-Froude-number  limit.  More 
precisely,  the  hull-surface  integral  can  be  approximated 
by  a  line  integral  around  the  mean  waterline  c. 
Furthermore,  this  integral  can  be  combined  with  the 
waterline  integral  in  equation  (26).  The  analysis  is 
presented  in  detail  in  Noblesse  (1986a)  and  is  briefly 
summarized  below 

The  mean  waterline  is  represented  by  the 


parametric  equations 

x  =  x(l(A)  and  y  -  y0(A),  (32a, b) 

where  the  parameter  A  varies  between  its  bow  and  stern 
values,  that  is 

AB  <  A  <  As  (32c) 

In  the  vicinity  of  the  mean  free  surface,  the  hull  surface 
is  represented  by  the  parametric  equations 

x  =  x0(A)  *  z  x [(A)  +  z*x2(A)  + .  .  .  ,  (33a) 

y  =  y0(A)  +  i  y,(A)  ♦  z2y,(A)  +  .  .  .  ,  <33b) 

where  l#  (  A  <  As  and  z  <  0.  (33c) 


The  velocity  potential  ^(A.z)  on  the  hull  surface  in  the 
vicinity  of  the  plane  z  =  0  likewise  is  expressed  in  the 
form 

♦  -  VA>  +  Z  ♦,(A)  +  zH2(A)v.  ..  .  (34) 

Differentiation  of  the  functions  xn( A) ,  y n( A),  ^(A)  with 
respect  to  the  parameter  A  is  denoted  by  the 
superscript  ;  thus,  we  have  x0  -  dx0(A)/dA 

By  using  the  foregoing  parametric  representations 
for  the  hull  surface  and  the  velocity  potential,  applying 


the  Laplace  method  for  approximating  the  hull-surface 
integral  in  equation  (26),  and  combining  the  resulting 
waterline  integral  with  that  already  present  in  equation 

(26),  we  may  obtain _ after  lengthy  algebraic 

transformations _ the  following  low-Froude-number 

approximation  for  the  function  K(t): 

K(t)  ~  v2^  J ^ (E0* a  t  +  E0_  a  _ )  dA  as  F  -  0,  (35) 

where  we  have 

q  =  1/p  =  1/(1  -i-t2)1'2,  (36) 

E0-  is  the  exponential  function 

Eg1  =  exp[-iv2p(x0±ty0)l,  (37) 

and  a±  is  the  amplitude  function 

at  =  uta1I  +  F2q2(uI)2a2t  +  0(F4),  (38) 

where  ut  is  defined  as 

ut  =  1/(1  -  iq(x, ity j )),  (39) 

and  the  functions  a,1  and  a,-  are  now  defined 

The  first-order  amplitude  function  a^  is  given  by 
a,1  =  yo'Aj/ll +t2)  +  2q(x0±ty0)(ut)2Bt40 

+  c±+o'  +  uDtf,/(l  +£2)  +  ip(y,V  -  y0>|>,  (40) 


where  we  have 

£  =  (y0 ' xi  -  *o'>i)/u-  (4I) 

with  u  defined  as 

u  =  f(x0')2  +  (y0')2l1'2,  (42) 

and  the  coefficients  A^,  B.,  C?  and  D,  are  defined  as 

A,  =  |(1  +pyn  'u)(l  -  py0'/u) +t2) 

+  dpy0'  u)[y,(x0'±ty0  )  u  +  t],  (43a » 

B.  =  q(y;*tx;)  +  i(y,x;  -  x,y,),  (43b) 

Ct  =  [<*t- p2x0y0'/u2) 


*  (py0‘  u)2t(x,x0'  +y,y0  )/u(l  +  e2)| 

♦  i(y,(x0  ±ty0')/u  +  c] 

|px0'  'u  -  (py0'/u)£(x,x0'  *  y,y0)  u(l  -  c2)].  (43c) 

D.  =  ((x0  ±ty0')/u)((l +£2)(y|±iqi)u1  +  i(py0  ulty,] 

-  (py0  /  u)£(py0' /u  —  if).  (43d) 

The  second-order  amplitude  function  a2*  is  given 

by 

a,4  y,  +  2uJ(+0m2't  •»  ^m,1 +  ^2m0*  iy0  r2*l 

*  bilUj^dglm^yj4  +  mg4/,1)  +  ♦)m01y21l 

I2(ut)\m01(y21)2.  (44) 


where  we  have 


Tn*  = 

(45) 

mn‘  =  H„  + 

(46) 

with 

Ho  =  £U- 

<47a) 

Hi  =  Xjy, '  -ViV  +  2<x2>o'  -*2*0 

(47b) 

H2  =  yi*2'  +  2(x25'|'  -V2xr> 

+  3(x3y0'  -y,xg'). 

(47c) 

In  the  particular  case  when  the  phase  of  the 
exponential  function  Eg1  is  stationary,  that  is  at  a  point 
(x0,y0,0)  where  we  have  x0'±ty0'  =  0,  the  first-order 
amplitude  function  a,1  takes  the  form 

a,1  =  ±  t|i>0  +<pyl+it)l0  +u+,l/(l±i£) 

+  ipiyiV  ■  y0 ' ♦  i >  if  V-%'  =  0  <48> 

Furthermore,  if  the  hull  surface  intersects  the  plane  z  =  0 
orthogonally  at  a  point  of  stationary  phase  we  have 

v  =  1  F:py0a\ax2 

if  x()'±tyg'  =  0  and  n,  =  0,  (49) 

and  the  amplitude  function  a,  then  is  or  order  E2. 

The  low-Froude-number  asymptotic  approximation 
to  the  far-field  wave-amplitude  function  K(t )  given  by  the 
waterline  integral  (35)  is  considerably  simpler  than  the 
exact  expression  (26),  which  involves  both  a  waterline 
integral  and  a  hull-surface  integral;  and  the  approximate 
expression  (35)  is  well  suited  for  efficient  numerical 
evaluation.  However,  expression  (35)  can  be  simplified  by 
applying  the  method  of  stationary  phase,  which  takes 
advantage  of  the  rapid  oscillations  of  the  exponential 
function  Fnr  given  by  equation  (A’)  in  ihe  low-Froude- 
number  limit  v  -*  »,  or  more  generally  tn  the  limn 
vp  This  stationary -phase  approximation  is  now 

obtained 

5  I  OW  I  ROI  DF  Nl  MBFR  STATIONARY  PHASE 
APPROXIMATION  TO  THE  FAR  I  II  I  I)  WAV!-- 
AMPl  IFI  DP  PL  NOTION 

The  method  ol  stationary  phase  indicates  that  the 
major  contributions  to  the  integral  ( A 5 )  in  the  limit  when 
the  exponential  (unctions  p(l'  are  rapidly  oscillating,  that 
is  if  v  —  oo  or  more  generally  vp  —  °°,  stem  from  points 
where  the  phases  of  these  exponential  functions  arc 
stationary,  that  is  Irom  points  where  we  have  xn'*tyn 
0,  and  from  the  end  points  Ag  and  As  ot  the  integration 


range,  that  is  from  the  bow  and  the  stern.  This 
stationary-phase  analysis  is  presented  in  detail  in 
Noblesse  (1986b)  and  only  its  re-  jits  are  given  here.  We 
have 

K(t)  ~  iq3(KB  -Ks  +  IKt)asF-  0,  (50) 

where  KB  and  Ks  correspond  to  the  contributions  of  the 
bow  and  the  stern,  respectively,  and  corresponds  to 
the  contribution  of  a  point  where  the  phase  of  the 
exponential  function  Eg1  is  stationary,  that  is  where  we 
have  y0' /( -  x0  )  =  ±l/t  or  dy0/dx0  =  *l/t;  the 
summation  in  the  asymptotic  approximation  (50)  is 
extended  over  all  such  points  of  stationary  phase  on  the 
mean  waterline  c.  The  expressions  for  the  contributions 
of  the  bow  and  stern  KB  s  and  of  the  points  of 
stationary  phase  K,  are  given  below. 

The  contribution  KB  s  of  the  bow  or  stern  may  be 
expressed  in  the  form 

KB  s  =  A„  sexp(  -  iv2pxg  s).  (51) 

where  xB  s  is  the  abscissa  of  the  bow  or  stern  and  the 
amplitude-function  AB  s  is  given  by 

Ab  s  c  2nxA,  *  iF2q(A;  ‘  +  A,  )  +  0(F4);  (52) 

in  this  expression,  the  functions  A!  and  A,!  are  defined 
by  the  equations 

(nx  +  iqn/)A)  -  nxnx/l  -  ipn,(nx  +  iqn,)/(l  -  n,2  -  p2nx2)) 
*!,♦,  •ipnxtx*J(l  n,2*ipnxn,)  (1  n?2  -  p2nx2)  (53) 

and 

9.  A,:  -  ( u , a ,  *  9^  )  *  iqtu,  )2a,',  (54) 

where  the  superscript  denotes  differentiation  with 
respect  to  the  parameter  A,  the  functions  u,,  a^  and  a-," 
are  defined  bv  equations  (49).  (40)  and  (44).  respectively, 
and  the  function  9.  is  defined  as 

8. 

Expression  (54)  for  the  second-order  amplitude- 
function  As*  is  a  complex  one.  However,  equation  (53) 
defines  the  first-order  amplitude  function  A,  explicitly  in 
terms  of  the  value  of  t.  p  (I  ♦  t2)1  2  and  q  =  1  p,  the 
geometrical  characteristics  of  the  hull  at  the  bow  or 
stern,  namely  the  unit  vector  t(t(,tx,0)  tangent  to  the 
mean  waterline  and  the  unit  vector  ri*(nx,nv,n;)  normal 
to  the  hull,  and  the  components  4,  and  of  the  velocity 
vector  in  the  directions  of  the  unit  vectors  Tand  n  *T 
tangent  to  the  hull.  Equation  (52)  shows  that  the  lirst 
order  approximation  to  the  amplitude  function  ABV 


given  by  2nxA,,  vanishes  if  n^-  0,  that  is  if  the 
waterline  has  a  cusp  at  the  bow  or  stern.  Equation  (S3) 
shows  that  the  first-order  approximation  to  the  function 
AB  s  also  vanishes  if  the  bow  or  stern  is  round,  since  we 
then  have  0^  =  0  and  ♦,  =  -  4y  =  0  by  symmetry.  It 
may  thus  be  seen  that  the  contribution  of  the  bow  or 
stern  to  the  far-field  wave-amplitude  function  K(t)  is  of 


order  F2  if  nx  =  0  or  ny  =  0  at  the  bow  or  stem, 
respectively;  that  is,  we  have 

kB,S  =  W2)  if  nx  =  0  or  ny  =  0  (56) 

at  the  bow  or  stern.  In  the  particular  case  when  the  hull 
surface  is  vertical  at  the  bow  or  stern  we  have  nz  =  0 
and  equation  (S3)  becomes 

A,  =  tx  -  +t  -  ipt„4z/(l  —  p^ty2)  if  nz  =  0,  (57) 

which  yields 

Ai  =  t,-+,  +  0(F2)  if  nz  =  0  (58) 


as  is  indicated  by  the  free-surface  boundary  condition  +z 

=  -?\x  =  o(f2). 

The  contribution  Kt  of  a  point  of  the  mean 
waterline  c  where  the  phase  of  the  exponential  function 
Eq1  is  stationary,  that  is  where  we  have 

dy0/dx0  =  Tl/t,  (59) 

may  be  expressed  in  the  form 
Kt  =  ±v(2nr)|/2A± 

exp  [siv^ypt,  -  x0ly)/ty2  +  icn/4],  (60) 

where  r  is  the  radius  of  curvature  of  c  at  the  point  of 
stationary  phase  (x^y,,).  i  is  equal  to  +  1  or  -  lif  the 
center  of  curvature  of  c  at  the  point  (Xp.yp)  is  upstream 
or  downstream  from  (x0,y()),  respectively,  and  the 
amplitude-function  \±  is  given  by 

At  =  (I  -nzJ)'/2(nzty  +  ♦,,)-  F2qA2*  +  0(F4);  (61) 

in  this  expression,  the  second-order  amplitude-function 
Aj1  is  defined  by  the  equation 

±2((x0)2  +  (y0')2ll/2A21  = 

+  u±a , 1  (5(9t )2/38t  "  -  0J4>1  /4<8  z ' '  )2 
+  2iq(u±)2a2I,  (62) 

where  the  superscript  '  denotes  differentiation  with 
respect  to  the  parameter  K  and  the  functions  ut,  a,1,  a2* 
and  8t  are  given  by  equations  (39),  (40),  (44)  and  (55), 
respectively. 


The  expression  for  the  second-order  amplitude- 
function  A2*  is  a  complex  one.  However,  equation  (60) 
and  the  first-order  approximation  to  the  amplitude 
function,  namely 

At  -  (1  —  nz2)l/2(n,ty  +  ♦„) ♦  Off1).  (63) 

provide  a  simple  explicit  expression  for  the  stationary- 
phase  contribution  K±  in  terms  of  the  geometric 
characteristics  of  the  hull  and  the  downward  tangential 
derivative  +d  of  the  potential  at  the  point  of  stationary 
phase.  In  the  particular  case  when  the  hull  surface  is 
vertical  at  the  point  of  stationary  phase  we  have  1^*0 
and  equation  (63)  becomes 

A±  =  -♦j  +  CXF2)  if  nz  =  0,  (64) 

which  yields 

A±  =  0(F2)  if  nz  =  0.  (65) 

Equation  (60)  then  shows  that  the  stationary-phase 
contribution  at  a  point  (x^yg)  of  c  where  the  hull  is 
vertical  is  of  order  F,  that  is  we  have 

K±  =  O(F)  if  nz  =  0.  (66) 

On  the  other  hand,  equations  (60)  and  (63)  show  that  we 
have 

K±  =  0(I/F)  if  nz  **  0.  (67) 

The  stationary-phase  contribution  at  a  point  where  the 
hull  has  flare  thus  is  dominant  in  the  zero-Froude- 
number  limit. 

The  summation  in  equation  (50)  is  extended  to  all 
the  points  of  the  mean  waterline  c  where  the  phase  of 
the  exponential  function  E0+  or  the  function  E0  is 
stationary,  that  is  the  points  where  the  slope  dy0/dxQ  of 
c  is  equal  to  -  l/t  or  +  1/t,  respectively.  The  number  of 
stationary  points,  and  their  position  along  the  waterline, 
depend  on  the  value  of  t  and  on  the  shape  of  c.  For 
instance,  for  the  simple  case  of  a  hull  with  waterline 
consisting  of  a  sharp-ended  parabolic  bow  region  1/4  < 
x  <  1/2  defined  by  the  equation  y  =  4bx(l  -2x),  where 
b  denotes  the  ship’s  beam/length  ratio,  a  straight  parallel 
midbody  region  -  1/4  <  x  <  1/4,  and  a  round-ended 
elliptic  stern  region  -  1/2  <  x  <  -  1/4  defined  by  the 
equation  y  =  b[-2x(l  +  2x))1/2,  there  is  one  point  of 
stationary  phase  in  the  stem  region  given  by  x  = 

-  1 1  +  1/(1  +4b2t2)1/2]/4,  so  that  we  have  -  1/2  <  x  < 

-  1/4  for  0  <  t  <  00  with  x  -*  -  1/2  as  t  -»  0  and  x  -» 

-  1/4  as  t  -»  oo,  and  one  point  of  stationary  phase  in  the 


bow  region  given  by  x  =  (1  +  l/4bt)/4  for  l/4b  <  t  < 

«°.  so  that  we  have  1/2  >  x  >  1/4  for  l/4b  <  t  <  00 
with  x  -•  1/2  as  t  -*  l/4b  and  x  -»  1/4  as  t  -»  We 
thus  have  one  point  of  stationary  phase  in  the  stern 
region  for  0  <  t  <  l/4b  and  two  points  of  stationary 
phase,  one  in  the  stern  region  and  one  in  the  bow  region, 
for  l/4b  A  t  <  <».  The  two  points  of  stationary  phase 
approach  the  shoulders  x  =  ±1/4,  where  dy/dx  =  0,  as 

t  -*  <*>. 

The  asymptotic  approximation  (50)  and  equations 
(51  )-(58)  and  (59)-(67)  defining  the  contributions  of  the 
bow  and  stern  and  of  the  stationary-phase  point(s)  on 
the  waterline,  respectively,  show  that  the  low-Froude- 
number  behavior  of  the  far-field  wave-amplitude 
function  is  strongly  influenced  by  the  shape  of  the  hull  in 
the  vicinity  of  the  waterline.  More  precisely,  for  a  value 
of  t  for  which  there  is  one  (or  more)  point  of  stationary 
phase  on  the  mean  waterline  where  the  hull  has  flare,  the 
contribution  of  this  stationary-phase  point  dominates  the 
contribution  of  the  bow  and  stern  and  is  of  order  1/F, 
that  is  we  have  K(t)  =  0(1/F)  On  the  other  hand,  for  a 
value  of  t  for  which  cither  there  corresponds  no  point  of 
stationary  phase  or  the  hull  has  no  flare  at  the  point(s) 
of  stationary  phase,  the  dominant  contribution  stems 
from  the  bow  and  stern,  and  it  is  of  order  1,  that  is  we 
then  have  K(t)  -  0(1).  However,  if  n^  =  0  or  ny  «  0 
at  both  the  bow  and  the  stern,  that  is  if  the  bow  and  the 
stern  are  either  cusped  or  round,  their  contribution  is 
<)<!-“)  and  the  contribution  of  the  stationary-phase 
pointls),  which  is  O(f)  if  there  is  no  flare  (as  is  assumed 
here),  is  dominant;  so  that  we  then  have  K(t)  =  O(F). 

I  or  a  ship  form  that  is  everywhere  wall  sided,  the 
contribution  of  the  bow  and  stern  is  dominant  for  all 
values  of  t,  and  we  have  K(t)  -  0(1)  for  0  <  I  <  ®.  On 
the  other  hand,  for  a  hull  form  that  has  flare  over  a 
portion  of  the  waterline  and  is  wall  sided  elsewhere,  the 
contribution  ot  the  bow  and  stern  is  dominant,  and  0(1), 
only  for  those  values  of  t  for  which  the  corresponding 
points  of  stationary  phase  fall  outside  the  range  of  flare; 
loi  the  range  ol  values  ot  t  for  which  the  corresponding 
points  of  stationary  phase  are  within  the  range  of  flare, 
the  contribution  of  these  stationary-phase  points  is 
dominant,  ot  order  1/F.  In  this  instance,  the  function 
kit)  is  0(1  I  )  for  a  range  of  values  of  t  (corresponding 
to  the  region  of  flare)  and  0(1)  for  other  values  to  t.  If 
the  region  of  flare  is  of  small  extent,  and  the  slope 


region,  the  range  of  values  of  t  for  which  K(t)  =  0(1/F) 
is  also  small.  The  far-field  wave-amplitude  function  K(t) 
can  then  exhibit  a  sharp  peak  for  some  value  of  t  in  the 
low-Froude-number  limit.  In  fact,  several  isolated  peaks 
of  the  function  K(t)  can  exist  if  the  hull  form  has  several 
distinct  regions  of  flare  within  which  the  slope  of  the 
waterline  varies  gradually. 

It  should  be  noted  that  the  result  K(t)  =  0(l/F) 
for  values  of  I  for  which  the  hull  has  flare  at  the 
corresponding  points  of  stationary  phase  does  not  imply 
that  the  corresponding  free-surface  elevation  becomes 
unbounded  in  the  zero-Froude-number  limit  F  =  0. 
Indeed,  the  asymptotic  approximation  (6),  where  we  have 
(«)  =  (X.Z)g/U2,  then  yields  Zg/U2  =  0[1/F(  ()l/2)  = 
0[l/(-X/L)1/2l  as  F  -  0  and  -  X/L  -  °°.  The  free- 
surface  elevation  Z  thus  is  of  order  (U2/g)/(  -  X/L)|/2  as 
F  -»  0  and  -X/L  *•  <*>,  and  equation  (15)  shows  that 
the  corresponding  wave  steepness  s  is  0|l/(  -  X  T.)1  2] 

For  values  of  t  for  which  K(t)  =  0(1)  it  is  seen  that 
Zg/U2  and  s  are  0(F/(  -  X/L)1  2)  as  F  —  0  and  -  X/L 

— *  OO 

It  should  also  be  noted  that  the  asymptotic 
approximation  (5 1 )-( 5 5 )  for  the  contribution  ol  the  bow 
and  stern  is  not  uniformly  valid  for  the  values  of  t  for 
which  the  bow  or  the  stern  is  a  point  of  stationary 
phase,  that  is  for  which  the  waterline  slope  dy0/dxn  at 
the  bow  or  the  stern  is  equal  to  I '  t  or  l/t, 
respectively.  Indeed,  it  may  be  shown  that  we  have 
1  -  n72-  p2nx2  =  0  at  a  point  of  stationary  phase,  so 
that  the  first-order  approximation  to  the  amplitude 
functions  AB  s  given  by  equation  (53)  becomes 
unbounded.  The  asymptotic  approximation  (60M62)  for 
the  contribution  of  a  point  ol  stationary  phase  likewise  is 
not  uniformly  valid  at  a  stationary'  phase  point  where  th 
waterline  has  an  inflexion  point  Indeed,  the  radius  of 
curvature  r  at  such  an  inflexion  point  is  infinite  and 
equation  (60)  yields  an  unbounded  contribution  k. 
Asymptotic  approximations  valid  for  these  special  cases 
may  easily  be  obtained  and  arc  given  in  Noblesse 
(1986b).  It  will  only  be  noted  here  that  the  far  held 
wave-amplitude  function  k(t)  for  a  hull  form  having 
flare  at  a  point  where  the  waterline  has  an  inflexion  may 
be  expected  to  exhibit  a  particularly  pronounced  peak  at 
the  value  of  t  corresponding  to  the  inflexion  point  since 
we  have  K(t)  =  0(1  T4  ’)  as  I  -*  0  for  this  particular 
value  of  t. 


dy  dx  of  the  waterline  does  not  vary  widely  within  that 


6.  CONCLUSION:  HULL  FORM  AND  KELVIN- 
WAKE  FEATURES 

The  classical  far-field  asymptotic  approximation  to 
the  Kelvin  wake,  obtained  in  section  2  by  applying  the 
method  of  stationary  phase,  and  the  low-Froude-number 
asymptotic  approximation  to  the  far-field  wave-amplitude 
function,  obtained  in  sections  4  and  5  by  successively 
using  the  Laplace  method  and  the  method  of  stationary 
phase,  provide  a  simple  analytical  relationship  between 
the  hull  shape,  on  one  hand,  and  the  waves  it  generates, 
on  the  other  hand.  This  explicit  relationship  between  the 
wavemaker  and  its  waves  is  summarized  below. 

The  far-field  asymptotic  approximation  to  the 
Kelvin  wake  (6)  shows  that  at  any  point  (i.a).  with  t  = 
Xg/U2  «  -  I  and  0  <  a  <  tan  “  '(2 “  3/2)  ~  19°28\  the 

wave  field  consists  in  two  plane  progressive  waves _ a 

transverse  wave  and  a  divergent  wave _ with 

wavelengths  A  _  and  A  t  propagating  at  angles  /3  _  and 
P from  the  track  of  the  ship,  respectively.  The 
wavelengths  AI  and  the  propagation  angles  px  depend  on 
the  angle  from  the  ship  track  a  alone,  that  is  X±  and  ft± 
are  independent  of  the  hull  shape  and  size,  as  is  well 
known.  Specifically,  the  functions  AI(o)  and  pjo)  are 
defined  by  equations  (10),  (11)  and  (12).  At  a  given 
downstream  distance  {,  the  amplitudes  of  these  waves, 
on  the  other  hand,  are  given  by  the  product  of  the 
functions  A  Jo),  defined  by  equation  (9),  and  the  far- 
field  wave-amplitude  function  K(t)  evaluated  at  the 
stationary  values  t±(o)  given  by  equation  (7).  The 
function  K(t)  depends  on  the  hull  shape  and  the  Froude 
number  in  a  fairly  complicated  manner  via  an  integral 
over  the  mean  wetted-hull  surface  and  an  integral  around 
the  mean  waterline,  as  is  indicated  by  equation  (26).  A 
low-Froude-number  asymptotic  approximation  to  these 
integrals  is  obtained  in  sections  4  and  5. 

The  analytical  approximation  (50)  shows  that  for  a 
given  value  of  t  corresponding  to  a  given  value  of  a,  as 
is  specified  by  equation  (7),  the  main  contributions  to  the 
function  K(t)  stem  from  several  particular  points  on  the 
mean  waterline.  These  are  the  bow  and  the  stern,  on  one 
hand,  and  (usually  but  not  always)  one  (or  several) 
point(s)  of  stationary  phase.  Indeed,  the  number  of  these 
points  of  stationary  phase  and  their  position  on  the 
waterline,  defined  by  the  condition 

|dy0/dx0|  =  1/t,  (68) 

depend  on  the  shape  of  the  waterline  and  the  value  of  t. 


Let  <fi  denote  the  angle  between  the  tangent  to  the 
mean  waterline  and  the  track  of  the  ship,  that  is  we  have 

tan  ip  *=  |dy0/dxg|  and  0  <  <p  <  n/2.  (69) 

Equations  (68),  (69),  (7)  and  (10)  then  show  that  the 
stationary  values  t  +  (o)  and  t.(o)  associated  with  a  given 
value  of  a  are  defined  by  the  equation 

tfijo)  =  tan  _  '{4o/[l±(l  -  8o2)l/2]},  (70) 

where  o  =  tana.  For  a  given  waterline  shape,  equation 
(70)  thus  defines  the  number  of  stationary  points  and 
their  position  on  the  waterline  corresponding  to  any 
given  angle  a  inside  the  Kelvin  wake.  In  particular, 
equation  (70)  yields 

0  <  <p+  <  tan-'(2,/2)  ~  54°44'  <  c p_  <90°  (71) 

and  *v  2o,  <p_  ~  n/2 -a  as  o  -»  0.  (72a,  b) 

Points  of  the  waterline  with  slope  between  0  and  54’44 
thus  contribute  mostly  to  the  system  of  divergent  waves 
while  waterline  slopes  between  54°44  and  90“  mostly 
contribute  to  the  transverse  waves. 

Equations  (10),  (II),  (12)  and  (70)  define  the  wave¬ 
lengths  At,  the  wave-propagation  angles  p±  and  the 
waterline-tangent  angles  ipx  corresponding  to  a  given 
angle  from  the  ship  track  a.  The  functions  A„(a),  pjo) 
and  q>±(a)  are  depicted  in  Figure  3a,  where  the  subscripts 
T  and  D  are  used,  instead  of  -  and  +  ,  to  refer  to  the 
transverse  and  divergent  waves,  respectively.  The 
foregoing  relationships  between  a  and  A.  ft.  <p  may  be 
used  for  determining  the  angle  from  the  ship  track  a,  the 
wavelength  A  and  the  wave-propagation  angle  P 
corresponding  to  a  given  waterline-tangent  angle  <p. 
Specifically,  we  may  obtain  the  remarkably  simple 
relations 


a  =  tan'  1  [tan<p/(2  +  tan2q>)], 
A/2n  =  sin2cp  and  P  =  n/2  -  i p. 


(74a, b) 


The  functions  a(q>),  A(<p)  and  P(<p)  are  depicted  in  Figure 
3b.  Alternatively,  the  foregoing  relationships  among  <p.  o. 
A  and  p  can  be  represented  in  the  form  of  Figures  3c  and 
3d,  which  depict  the  functions  <p(P).  o(p).  A(/))  and  ip(A), 
o(A).  (3(A),  respectively.  These  equivalent  graphical 
representations  show  that  we  have 

0  <  tpD  <  tan  -  *<2l/2)  a#  54°44  <  <pT  <  90°,  (75a) 

0  <  AD/2n  <  2/3  <  AT/2n  <  I.  (75b) 

90°  >PD  >  tan-'(2-,/2)  a,  35°  16'  >  fty  >  0,  (75c) 
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Fig.  3a  -  The  Waterline-Tangent  Angle  cp,  the  Wave- 
Propagation  Angle  p  and  the  Wavelength  A  as  Functions 
of  the  Angle  From  the  Ship  Track  a;  the  Subscripts  T 
and  D  Refer  to  the  Transverse  and  Divergent  Waves, 
Respectively. 
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fig  3c  The  Waterline-Tangent  Angle  <p,  the  Angle 
from  the  Ship  Track  a  and  the  Wavelength  A  as 
functions  of  the  Wave-Propagation  Angle  ft. 
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Fig.  3b  -  The  Angle  From  the  Ship  Track  a,  the  Wave- 
Propagation  Angle  p  and  the  Wavelength  A  as  Functions 
of  the  Waterline-Tangent  Angle  <p. 
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Fig.  3d  The  Waterline-Tangent  Angle  ip,  the  Angle 
f  rom  the  Ship  Track  a  and  the  Wave-Propagation  Angle 
P  as  functions  of  the  Wavelength  A 


where  the  subscripts  D  and  T  refer  to  the  divergent  and 
transverse  waves,  respectively. 

The  first  two  terms  in  the  low-Froudc-number 
asymptotic  expansions  for  the  contributions  KH  s  of  the 
bow  and  stern  and  the  contributions  K.  of  the  points  of 
stationary  phase  in  equation  (50)  are  given  by  equations 
(5l)-(55)  and  <59)-(62),  respectively  The  second-order 
terms  in  these  asymptotic  expansions  are  defined  by 
complex  expressions  However,  the  first  order  terms 
provide  simple  approximations  defined  explicitly  in  terms 
of  the  geometrical  characteristics  of  the  hull  and  the 
velocity  components  in  the  tangential  directions  T and 
nVTto  the  hull.  In  particular,  the  low-Froude-number 
asymptotic  expansions  given  in  section  5  show  that  the 
contributions  K0  and  Ks  of  the  bow  and  stern  are  0(1) 
except  if  the  bow  or  stern  is  cusped  or  round,  in  which 


case  we  have  KH  s  =  Of f  ‘ )  The  contribution  of  a  given 
point  of  stationary  phase  is  0(1 /F).  and  thus  is 
dominant,  if  the  hull  has  flare  at  this  point;  otherwise, 
that  is  if  the  hull  is  wall  sided  at  the  point  of  stationary 
phase,  its  contribution  is  O(F) 

Thus,  for  a  ship  form  that  is  everywhere  wall 
sided,  the  contributions  of  the  bow  and  stern  are 
dominant  (assuming  that  they  are  not  both  either  round 
or  cusped)  for  all  values  of  the  angle  o  from  the  ship 
track,  that  is  everywhere  in  the  Tar-field  Kelvin  wake.  On 
the  other  hand,  for  a  hull  that  has  flare  over  a  portion 
of  the  waterline  and  is  wall  sided  elsewhere,  the 
contributions  of  the  bow  and  stern  are  dominant,  and 
0(1),  only  for  those  angles  a  in  the  Kelvin  wake  for 
which  the  corresponding  points  of  stationary  phase  on 
the  waterline  fall  outside  the  region  of  flare;  for  the 


range  of  values  of  a  for  which  the  corresponding  points 
of  stationary  phase  are  within  the  range  of  flare,  the 
contribution  of  these  stationary-phase  points  is 
dominant,  of  order  l/F.  In  this  instance,  the  amplitude 
of  the  waves  in  the  far-field  Kelvin  wake  is  of  order 
(U2/g)/(  -  X/L)1'2,  as  F  -*  0  and  X/L  -*  for  the 
range  of  values  of  a  corresponding  to  the  region  of  flare 
and  of  order  F(U2  g)/(  -  X/L)1  2  for  values  of  a  outside 
this  range.  If  the  region  of  flare  is  of  small  extent,  and 
the  waterline-tangent  angle  does  not  vary  widely  within 
that  region,  the  corresponding  range  of  values  of  the 
angle  o  where  the  wave  amplitude  is  an  order  of 
magnitude  larger  than  elsewhere  is  also  small,  and  thus 
appears  as  a  peak  for  sufficiently  small  values  of  the 
Froude  number.  This  peak  is  particularly  pronounced  for 
a  hull  with  a  small  region  of  flare  in  the  vicinity  of  an 
inflexion  point  of  the  waterline. 

The  low-Froude-number  asymptotic  analysis  of  the 
Neumann-kelvm  theory  presented  in  this  study  thus 
shows  that  the  characteristics  of  the  far-field  Kelvin  wake 
strongly  depend  on  the  shape  of  the  ship  hull,  notably 
the  presence  of  flare  and  the  shape  of  the  waterline  at 
the  bow  and  stern  This  analysis  also  predicts  that  the 
nondinicnsional  wave-resistance  coefficient  R/fiU2L“, 
where  l  and  I  are  the  speed  and  the  length  of  the  ship 
and  y  is  the  density  of  water,  which  is  given  by  the 
Havelock  integral 

nR  yl  '1  2  ,4/0  K(t)  2(l  -M2)  1  2  dt,  (76) 

is  0(1  ’)  lot  a  ship  form  with  a  region  of  flare,  Q(F4)  for 
a  ship  lorm  that  is  wall  sided  everywhere  and  has  either 
a  bow  or  a  stern  (or  both)  that  is  neither  cusped  nor 
loiind.  and  (>()*)  for  a  wall-sided  ship  form  with  both 
bow  and  stern  that  arc  cither  cusped  or  round 

the  low -t  roude-numher  asymptotic  approximation 
(50)  also  shows  that  we  have 

Km  0(1  t’)  as  t  —  °°  (77) 

This  result  is  actually  valid  for  any  value  of  the  Froude 
number,  indeed,  the  asymptotic  approximation  (50)  is 
valid  not  only  in  the  low-Froude-number  limit  F  0  but 
more  generally  in  the  limit  vp  (I  -ft2)1  2  F  -*  <*>,  that 
is  as  F  —  0  or/and  as  t  —  <*>,  as  may  be  seen  from  tne 
exponential  functions  exp  (v2p2i)  and  Et  - 
exp  |  iv2p(x-tty)|  in  equation  (26) 

Equations  (15),  (18a),  (21b)  and  (77)  then  yield 

(  4)1  2s  (f.u)  ~  t  #  1  2  ~  I  o1  2  as  a  -*  0  (78) 


for  the  steepness  of  the  short  divergent  waves  in  the 
vicinity  of  the  track  of  the  ship.  By  using  equations  (4) 
and  (5),  which  yield  a  ~  iy/(  -  i)  as  a  -*  0,  in  equation 
(78)  we  may  obtain 

s  t  (t.rj)  ~  1/rj1  2  as  r)  -»  0.  (79) 

Equation  (79)  thus  shows  that  the  lines  along  which  the 
steepness  of  the  short  divergent  waves  in  the  far-field 
Kelvin  wake  takes  given  large  values,  say  s(  =  1/7  and 
1/15,  are  parallel  to  the  ship  (rack,  as  was  found  in 
Figure  21  of  Barnel)  and  Noblesse  (1986)  by  using  the 
Michel!  thin-ship  approximation  for  a  simple  ship  form. 
The  Neumann-Kelvin  theory  therefore  predicts  that  the 
far-field  Kelvin  wake  contains  three  distinct  regions:  (i)  a 
narrow  constant-width  inner  region  bordering  the  track 
of  the  ship  where  no  divergent  gravity  waves  can  exist, 

(n)  an  outer  region  where  the  usual  transverse  and 
divergent  waves  are  present,  and  (hi)  an  intermediate 
region  at  the  boundary  between  the  inner  and  outer 
regions  where  short  steep  divergent  waves  can  be  found. 
In  reality,  surface  tension  must  evidently  be  taken  into 
account  in  the  vicinity  of  the  track  of  the  ship. 
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